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Abstract. The Poisson-Boltzmann mean- field description of ionic solutions has been 
successfully used in predicting charge distributions and interactions between charged 
macromolecules. While the electrostatic model of charged fluids, on which the Poisson- 
Boltzmann description rests, and its statistical mechanical consequences have been 
scrutinized in great detail, much less is understood about its probable shortcomings 
when dealing with various aspects of real physical, chemical and biological systems. 
These shortcomings are not only a consequence of the limitations of the mean- 
field approximation per se, but perhaps are primarily due to the fact that the 
purely Coulombic model Hamiltonian does not take into account various additional 
interactions that arc not electrostatic in their origin. We explore several possible 
non-electrostatic contributions to the free energy of ions in confined aqueous solutions 
and investigate their ramifications and consequences on ionic profiles and interactions 
between charged surfaces and macromolecules. 
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1. Introduction 

The traditional approach to ions in solution has been the mean-field Poisson-Boltzmann 
(PB) formalism. This approach adequately captures the main features of electrostatic 
interactions at weak surface charges, low ion valency, and high temperature [HE]- It 
stems from a Coulomhic model Hamiltonian that includes only purely electrostatic 
interactions between different charged species. The limitations of the PB approach, 
which rests on a collective and continuous description of statistical charge distributions 
become particularly important in highly-charged systems, where counterion-mediated 
interactions between charged bodies cannot be described by the mean-field approach 
that completely neglects ion correlations and charge fiuctuations [3]. These mean-field 
limitations have been successfully bypassed and have led to more refined descriptions 
that capture some of the important non- mean- field aspects of Coulomb fluids [1]. 

Beside electrostatic interactions that are universal, omnipresent non-electrostatic 
interactions are specific and dependent on the nature of the ionic species, solvent and 
confining interfaces. Because it is difficult, perhaps impossible, to devise a universal 
theory accounting for all non-electrostatic effects, such additional interactions need to 
be treated separately, describing specific electrolyte features that go beyond regular PB 
theory. These effects depend on ionic chemical nature, size, charge, polarizability and 
solvation (preferential ion-solvent interaction ^). 

In this review we describe how adding non-electrostatic terms in the system free 
energy yields modified PB equations and ionic density profiles and differentiate between 
the way different ionic species interact with macromolecules. While it has long been 
recognized that ions can have such different effects on macromolecular interactions [8], 
only recently have such ionic features been accounted for within electrostatic mean- 
field theory. Examples include ion effects on surface tension [9] and precipitation of 
proteins from solutions [lO]. The latter led to the so-called Hofmeister ranking of 
different ions according to their surface activity. Adsorption of ions and/or surfactants to 
charged interfaces also contains substantial contributions from non-electrostatic degrees 
of freedom, as was shown within the Ninham-Parsegian theory of charge regulation at 
surfaces [TT] . 

The outline of this paper is as follow. After reviewing the standard PB theory in 
section 2, we present in section 3 how addition of steric effects result in saturation of 
ionic profiles close to charged interfaces. In section 4 we show how non-electrostatic 
interactions between ions and charged membranes can cause a phase transition between 
two lamellar systems of different periodicity. Solvation effects are the topic discussed 
in section 5, where local variation of the dielectric function in solvent mixtures and 
ion-solvent interactions lead to changes in ionic and solvent profiles close to charged 
interfaces. Finally, in section 6 possible polarization effects of ions in solution is added 
to their ionic character, again resulting in different behavior close to charge interfaces. 



Ion-specific interactions in aqueous solutions 



3 



a 



a 




z=-D/2 



z=0 



z=+D/2 



Figure 1. Schematic illustration of the model system. The two plates residing at 
z = ±D/2 are charged with surface charge density cr. The electrolyte ions are denoted by 
© and 0. Their densities are n+{z) and 



2. The Poisson-Boltzmann model: summary and main results 

The PB theory is a useful starting point for many theoretical ramifications because 
it relies on a simple and analytically tractable model that can be easily extended 
and amended. Although the theory has its own well-understood limitations, it yields 
meaningful results in good agreement with many experiments [Tj. 

The model can be cast in many geometries but we shall focus on the simplest planar 
system, as depicted in figure [H Two charged planar plates (of infinite extent), located at 
z = ±D/2, are immersed in an electrolyte bath. Each plate carries a charge density of a 
(per unit area) and their inter-surface separation D is a tunable system parameter. The 
ionic solution is denoted clS Zj^ '. Z— and is composed of a solvent of dielectric constant e 
and two types of ions: cations of valency z+ and local density (number per unit volume) 
n+{z), and anions of valency 2_ with local density n-{z). The finite system of thickness 
D is in contact with an electrolyte bath of bulk densities obeying the electroneutrality 
condition: = Z-u'L- Throughout the paper, we will use the convention that Z- 

(denotes the valency of the anions) is taken as a positive number, hence, their respective 
charge is written as —ez_. 

The free energy within the PB model FpB can be derived in numerous ways as can 
be found in the literature [Tj. With our notation Fp-Q is a function of the densities Uj- 
and the local electrostatic potential ip: 
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- J d^r[n+n+ + fi^n^] (1) 

The above free energy is a functional of three independent fields: n± and ip- The first 
integral in ([1]) is the electrostatic energy, while the second represents the ideal mixing 
entropy of a dilute solution of the ± ions. Note that a microscopic length scale a was 
introduced in the entropy terms above. Only one such microscopic length scale will be 
used throughout this paper, defining a reference density associated with close-packing 
no = l/a^. 

The last integral in ([1]) is written in terms of the chemical potentials fi± of the ± 
ions. Alternatively, the chemical potential can be regarded as a Lagrange multiplier, 
setting the bulk densities to be = exp{fi±/ k-sT) / . 

The thermodynamic equilibrium state is given by minimizing the above free energy 
functional. Taking the variation of the free energy ([1]) with respect to n± yields the 
Boltzmann distribution of the ions in the presence of the local potential ip: 

™ : ± ez±^ + k^T ln(n±a=') - /i± = (2) 



wherefrom 



5n4 



n± = n_|_ exp {^eZj-ip/k-BT) (3) 



Similarly, taking the variation with respect to the potential yields the Poisson equation 
connecting with the ± densities: 

— r— = — V^V + e-2+n+ — ez^n^ = (4) 

dip 4:71 

Combining the Boltzmann distribution with the Poisson equation yields the familiar 
Poisson-Boltzmann equation: 

= [z+nle-''+^/''^^ - z_n'_e''-^/''^^] (5) 

that serves as a starting point for various extensions presented in the sections to follow. 

In addition to the volume contribution ([1]) of the free energy, we need to include 
a surface electrostatic energy term Fs, which couples the surface charge density a with 
the surface value of the potential ip{z= ± D/2) = ips- This surface term has the form 

d^r aips (6) 



A 

and takes into account the fact that we work in an ensemble where the surface charge 
density is fixed. Variation of FpB + Fg with respect to ips yields the well-known 
electrostatic boundary condition 

— (FpB + F,) = ^h-V^ + a = (7) 

dips 4:71 

wherefrom 

471 

n . = — -a (8) 
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where n is a unit vector normal to the surface. This boundary condition can also be 
interpreted as the electroneutrality condition since the amount of mobile charge should 
exactly compensate the surface charge. 

The above PB equation is valid in any geometry. But if we go back to the planar 
case with two parallel plates, it is straightforward to calculate the local pressure Ppb 
at any point z between the plates. In chemical equilibrium PpB should be a constant 
throughout the system. Namely, PpB is independent of the position z and can be 
calculated by taking the proper variation of the free energy with respect to the inter- 
plate spacing D: 

PpB = -^{^T + ksT (n+ + n_) (9) 

OTT 

The pressure is composed of two terms. The first is the electrostatic pressure stemming 
from the Maxwell stress tensor. This term is negative, meaning an attractive force 
contribution acting between the plates. The second term originates from the ideal 
entropy of mixing of the ions and is positive. This term is similar to an "ideal gas" van 
't Hoff osmotic pressure of the ± species. 

A related relation can be obtained from for one charged surface. Comparing 
the pressure Ppb calculated at contact with the charge surface, with the distal pressure 
{z — > oo) results in the so-called Graham equation used in colloid and interfacial science 

m 



2 £kBT , 



27r 



< + n'_ - < - ni) (10) 



where are the values of the counterion and co-ion densities calculated at the surface. 

Note that the osmotic pressure as measured in experiments can be written as the 
difference between the local pressure and the electrolyte bath pressure: 

n = PpB - A;bT « + ni) (11) 

We now mention two special cases separately: the counterion only case and the 
linearized Debye-Hiickel theory. 

2.1. The counterion only case 

In the counterion only case no salt is added to the solution and there are just enough 
counterions to balance the surface charge. For this case we have chosen arbitrarily the 
sign of the surface charge to be negative, a < 0. By setting n+ = 0, the free energy is 
written only in terms of the anion density in solution, n = n_ and Z- = z 



PB 



d^r 



{VipY " Gzmp + kBT{n\rL n) 

Stt uq 



(12) 



The local pressure PpB also contains only osmotic pressure of one type of ions (n_ 
apart from the Maxwell stress term. 
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2.2. The linearized Debye-Hiickel theory 

When the surface charges and potentials are small, eips ksT {ips ^ 25 mV at room 
temperature), the PB equation can be linearized and matches the Debye-Hiickel theory. 
For the simple and symmetric 1:1 monovalent electrolytes the PB equation reduces to: 

V^^ = Ad^^ (13) 

where the dimensionless potential is defined as = eip/k^T and the Debye screening 
length is Ad = ^/ek^T/Sire'^nh. The main importance of the Debye length Ad is to 
indicate a typical length for the exponential decay of the potential around charged 
objects and boundaries: "^{z) ~ exp(— ^/Ad)- The Ad length is about 3A for = IM 
of NaCl and about 1/im for pure water. It is convenient to express the Debye length 
in term of the Bjerrum length /b = e'^/ek-QT as Ad = ^/l/Snl^ni,. For water at room 
temperature Ib — 7A. 

Returning now to the general non-linear PB model, its free energy ([T]) can also be 
expressed in terms of the dimensionless potential \E'. For a 1:1 symmetric electrolyte we 
write it as: 

1 



FpB/kBT = J dV 



-(V^)^ + {n+ -n_)^ 



Stt/b 

-|- y d^r [n+ ln(n+/nb) + n_ ln(n_/nb) 

- 72+ - n_ + 2nb] (14) 

where the bulk contribution to the free energy is subtracted in the above expression. 

In the next sections we will elaborate on several extensions and modifications of 
the PB treatment. Results in terms of ion profiles and inter-surface pressure will be 
presented and compared to the bare PB results. 



3. Steric effects: finite ion size 



At sufficiently high ionic densities steric effects prevent ions from accumulating at 
charged interfaces to the extent predicted by the standard PB theory. This effect has 
been noted already in the work of Eigen [13] , elaborated later in [H] and developed into 
a final form by Borukhov et al. [15] and more recently in [16]. Steric constraints lead 
to saturation of ion density near the interface and, thus, increase their concentration in 
the rest of the interfacial region. This follows quite generally from the energy-entropy 
competition: the gain from the electrostatic energy is counteracted by the entropic 
penalty associated with ion packing. Beyond the mean-field {e.g., integral equation 
closure approximations), the correlations in local molecular packing clearly lead to ion 
layering and non-monotonic interactions between interfaces [T71 [18] . 

Steric effects lead to a modified ionic entropy that in turn gives rise to a modified 
Poisson-Boltzmann (MPB) equation, governing the distribution of ions in the vicinity of 
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charged interfaces. The main features of the steric effect can be derived from a lattice- 
gas model introduced next [T9] . Other possible approaches that were considered include 
the Stern layer modification of the PB approach, extensive MC simulations or numeric 
solutions of the integral closures relations pOj . 

We start with a free energy where the entropy of mixing is taken in its exact lattice- 
gas form, without taking the dilute solution limit. Instead of point-like particles, the 
co-ions and counter-ions are now modeled as finite-size particles having the same radius 
a. The free energy for a z+:z- electrolyte is now a modification of FpB dl]): 



*MPB 



61T 



+ ksT j dV [n+ ln(a^n+) + n_ ln(a^n_)] 

d^r (1 — a^n^ — a^nJ) ln(l — a^n^ — a^nJ) 



+ 



— / d^r [/i+n+ + 



(15) 

Taking the variation with respect to the three fields: n± and ip, yields the MPB 
equilibrium equations. We give them below for two cases: (i) symmetric electrolytes 
z = = z^] and, (ii) l:z asymmetric ones. 

In the former case the MPB equation is written as: 

STcezrih sinh.{e zip / k-gT) 
e 1 — h + hcosh.{ezip /k-oT) 
while the local ion densities are 

eyi]i{^ezip / k-QT) 



(16) 



n + n cosh.{ezip / k^T) 



:i7) 



where n = a^ln^ + n^) is the bulk volume fraction of the ions. Clearly, the ion densities 
saturate for large values of the electrostatic potential, preventing them from reaching 
unphysical values that can be obtained in the standard PB theory. 

In the latter case of l:z electrolytes, the MPB is obtained in the form 



„2 , 47r , 

y ip = (en+ — ezn. 

with the corresponding local densities 



71- 



zube 



n 



Ube 



ze^/kBT 



(19a) 



As discussed above, for large electrostatic potentials the densities saturate at finite 
values dependent on n, z and Ub- 

Figure [2] shows the ion density profile close to a single charged interface with fixed 
surface charge density a, in contact with an electrolyte bath of ionic density Ub and for 




Figure 2. Counterions density profile calculated for the MPB model. The solid line 
represents the standard PB model, while the dotted and dashed lines represent the MPB 
results with molecular sizes of 7.5A and lOA, respectively. In (a) the salt is symmetric 1:1, 
while in (b) it is asymmetric 1:4. Other parameters are: g = — 1/50 A^^, = O.IM , and 
e = 80. 



several ionic sizes a. As is clear from figure [2] the steric constraints limit the highest 
possible density in the vicinity of the charged surface and, thus, extend the electrostatic 
double layer further into the bulk, if compared to the standard PB theory. The extent 
of the double layer depends crucially on the hardcore radius of the ions, a. Furthermore, 
the valency of the counterions also affects the width of the saturated layer, as clearly 
demonstrated by the comparison between figure [2^ and figure Wp- 

The local pressure Pmpb can be calculated with analogy to Prb of ^ 

1.12 



^MPB = - + 



m 



a3 



1 + ^ ^^_e-e.+ ^/fcBT ^ ^^^ez.^/k^T^ 

+ 1 — n 



(20) 



and can be cast into the form of the contact theorem that relates the value of the 
equilibrium osmotic pressure to the values of the surface potential and its first derivative, 
yielding a generalization of the standard Graham equation used in the colloid science 



2 eksT 1 1 - 2a-^nb 
(7 ~ -m 21) 

This analytical expression is valid in the limit that the co-ions have a negligible 
concentration at the surface, ni 0. It is instructive to find that for a — >■ 0, expanding 
the logarithm in (!2T|) . we recover the Graham equation for the standard PB model (fTOj) . 
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4. PB and charge regulation in lamellar systems 
4.I. The model free energy and osmotic pressure 

In the previous sections the surface free energy Fg ([6]) was taken in its simplest form 
assuming a homogeneous surface charge, i.e. in the form of a surface electrostatic 
energy. We consider now a generalized form of Fg where lateral mixing of charged 
species is allowed within the surface and is described on the level of regular solution 
theory [I9]. 

Experimentally observed lamellar-lamellar phase transitions in charged surfactant 
systems |2T] provides an example where non-electrostatic, ionic-specific interactions 
appear to play a fundamental role. The non-electrostatic interactions are limited 
to charged amphiphilic surfaces confining the ionic solution. As evidenced in NMR 
experiments, ions not only associate differently with the amphiphile-water interface, 
but their binding may also restructure the interface they contact [22]. Computer 
simulations also indicate that the restructuring of the amphiphilic headgroup region 
should be strongly influenced by the size of the counter- ion [23]. Such conformational 
changes at the interface are possible sources of non-ideal amphiphile mixing, because 
non-electrostatic ion binding at the interface may effectively create two incompatible 
types of amphiphiles: ion-bound and ion-detached. 

We proposed a model [21] based on an extension of the Poisson-Boltzmann theory 
to explain the first-order liquid-liquid (L^ — » Lq,/) phase transition observed in osmotic 
pressure measurements of certain charged lamellae-forming amphiphiles [21]. Our 
starting point is the same as depicted in figure [H The free energy of the confined 
ions has several contributions. The volume free energy, FpB, is taken to be the same as 
the PB expression for the counterion only case f|T2l) . Because all counterions in solution 
originate from surfactant molecules, their integrated concentration (per unit area) must 
be equal in magnitude and opposite in sign to the surface charge density 

I'D/2 

2a = e n{z) dz (22) 

J ~D/2 

This is also the electroneutrality condition and can be translated via Gauss' law into the 
electrostatic boundary condition (in Gaussian units): ip'{D/2) = ip'^ = Ai^a /e, linking 
the surface electric field ip'^ with the surface charge density a. 

The second part of the total free energy comes from the surface free energy, Fg, 
of the amphiphiles residing on the planar bilayers. Here, we deviate from the Fs 
expression in ([6]) because we allow the surfactants on the interface to partially dissociate 
their counterion in the spirit of the Ninham-Parsegian theory of charge regulation 
The surface free energy Fg has electrostatic and non-electrostatic parts as well as a 
lateral mixing entropy contribution. Expressed in terms of the dimensionless surface 
area fraction r/s = a^a/e of charged surfactants and dimensionless surface potential 
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\l/s = eips/k^T, the free energy Fg is: 



S 9 



W I 2 



A 



+ r^^lnr/s + (1 - r/s) ln(l - r/^ 



(23) 



The first term couples the surface charge and surface potential as in ([6]), while the 
additional terms are the enthalpy and entropy of a two-component liquid mixture: 
charged surfactant with area fraction rjs and neutralized, ion-bound surfactants with 
area fraction 1 — ?7s. The dimensionless parameters and Xs are phenomeno logical, and 
denote respectively the counterion-surfactant and the surfactant-surfactant interactions 
at the surface. Here, < means that there is an added non-electrostatic attraction 
(favorable adsorption free energy) between counterions and the surface; the more 
counterions are associated at the surface, the smaller the amount of remaining charged 
surfactant. A positive Xs parameter represents the tendency of surfactants on the surface 
to phase separate into domains of neutral and charged surfactants. 

The total free energy Ftot is written as a functional of the variables '^{z)^ n{z), and 
a function of rjs, and includes the conservation condition, (l22l) . via a Lagrange multiplier, 
A.: 



Ftot[^,n;r]s] = Fpb + - A, 



r.D/2 

2 ■ 



rjs — a / n{z) dz 



(24) 



Next, we minimize Ftot with respect to the surface variable r^^, and the two continuous 
fields n{z), "^{z): dFtot/drjs = 5Ftot/5n{z) = 6Ftot/S'^{z) = 0, corresponding to three 
coupled Euler-Lagrange (EL) equations. The first one connects the surface charge 
density rjs with the surface potential "ifs 

- exp (As + as + XsVs - ^.) (25a) 



The second one is simply the Boltzmann distribution for the spatially-dependent ion 
density 

a^n{z) = exp (-A, + ^(2)) (256) 

and the last one is the standard Poisson equation 

^"(z) = AnlBuiz) (25c) 

In addition, the variation with respect to \E's gives the usual electrostatic boundary 
condition of the form 

nD/2) = K = ^ (26) 

The Lagrange multiplier, A^, acts as a chemical potential with the important difference 
that it is not related to any bulk reservoir, but rather to the concentration at the 
midplane, n(0). 

The non-electrostatic, ion-specific surface interactions govern the surface charged 
surfactant area fraction rjs, and has the form of a Langmuir-Frumkin-Davis adsorption 
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Figure 3. (a) The osmotic pressure 11 in units of 10^ Pascals; and, (b) the area fraction 
rjs — a^a/e of surface charges, as function of inter-lamellar spacing D for = —6, 
Xs — 12 and a — Sk. The Maxwell construction gives a coexistence between a phase with 
D ~ 39A and low 77s < 0.1, and another with D ~ 64A and ijs ^ 1- In (b) the two 
coexisting phases are denoted by squares and the dotted-dashed line shows the tie-line in 
the coexisting region. 
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Figure 4. Fit to the experimental osmotic pressure isotherm Il{D) of Ref. [21] on a 
log-log scale. The diamonds and squares are the data points for DDABr and DDACI, 
respectively, reproduced from Ref. [21]. The solid line is the best fit of the model to the 
phase transition seen for DDABr with — —7.4, Xs — 14.75 and a — 8 A. The fit also 
includes a hydration contribution (parameters for the form: Ilhyd — Ho exp(— Z3/Ahyd), 
with typical values: Ho = 2.37-10* Pa and Ahyd — 1-51 A). This contribution is particularly 
important at the low D region of the DDABr isotherm. A small amount of salt is added 
in the fits in the experiment (n^ = 0.5 mM). The dashed line is the fit to the DDACI (no 
transition) and all parameters are the same as for the solid line, except Ug — —3.4. 




Ion-specific interactions in aqueous solutions 



12 



isotherm Combining (125 6|) and fl25cp together is equivalent to the PB equation. 

Their solution, together with the adsorption isotherm (125 ap , completely determines the 
counterion density profile n{z), the mean electrostatic potential "^{z), and yields the 
osmotic pressure 11. We will apply this formalism to a specific example of lamellar- 
lamellar phase transitions next [2^ . 

4-2. The lamellar-lamellar transition in amphiphilic systems 

The typical isotherm Il{D) shown in figure [3^ exhibits a first-order phase transition 
from one free energy branch at large inter-lamellar separation D to another at smaller 
D — with a coexistence region in between. For given values of and Xs (chosen in 
the figure to be as = —6 and Xs = 12), and for large enough D (D > 64A), most 
counterions are dissociated from surfaces, rjg < 1, and the osmotic pressure follows the 
standard PB theory for (almost) fully dissociated surfactants. For smaller values of D 
{D < 39A), most counterions bind to the surface rjs < 0.1 and the isotherm follows 
another branch, characterized by a much smaller surface charge of only about 10% of 
the fully dissociated value. In intermediate D range (39A < D < 64A), the system is 
in a two-phase coexistence, the osmotic pressure 11 has a plateau and rj^ changes from 
one branch to the second (figure [3b). 

Our model is motivated by experiments [21] on the surfactant homolog series: 
DDACl, DDABr and DDAI The main experimental observation is reproduced in 
figure 4 along with our model fittings. When Cl~ serves as the counterion, as in 
DDACl, the osmotic pressure isotherm n(D) follows the usual PB result. When Br~ is 
the counterion, as in DDABr, one clearly observes a lamellar-lamellar phase transition 
from large inter-lamellar spacing of 60Ato small inter-lamellar spacings of about lOA. 
In addition, for the largest counterion, I~, as in DDAI, the lamellar stack cannot be 
swollen to the large D values branch . We can fit the experimental data by assigning 
different as and Xs values to the three homolog surfactants as can be seen in figure 4. 
Qualitatively, the different lamellar behavior can be understood in the following way. 
The Cl~ counterion is always dissociated from the DDA+ surfactant resulting in a PB- 
like behavior, and a continuous n(D) isotherm. For the Br~ counterion, the dissociation 
is partial as is explained above (figure 3) leading to a first-order transition in the isotherm 
and coexistence between thin and thick lamellar phases. Finally, for the DDAI, the I~ 
ion stays associated with the DDA"*" surfactant and there is no repulsive interaction to 
stabilize the swelling of the stack for any 11 value. More details about our model and 
the fit can be found in plj. 

Non-electrostatic interactions between counterion-associated and dissociated 
surfactants can be responsible for an inplane transition, which, in turn, is coupled to 
the bulk transition in the interaction osmotic pressure as can be clearly seen in figure [31 
This proposed ion-specific interactions are represented in our model by Xs and as- While 

I DDA stands for dodecyldimethylammonium and CI, Br and I correspond to chloride, bromide and 
iodine, respectively. 
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at present direct experimental verification and estimates for tlie proper Xs values are 
lacking, the conformational changes induced by the adsorbing ion, together with van der 
Waals interaction between adsorbed ions can lead to significant demixing. Furthermore, 
because larger ions are expected to perturb the surfactant-water interface to a larger 
extent, it is reasonable to expect that the value of Xs will scale roughly with the strength 
of surface- ion interactions, a^. 

We note that the Xs values needed to observe a phase transition, typically ~ 10 
(in units of k^T), are quite high These high values are needed to overcome the 
electrostatic repulsion between like-charged amphiphiles, leading to segregation. The 
source of this demixing energy, as codified by Xs, could be associated with mismatch of 
the hydrocarbon regions as well as headgroup-headgroup interactions, such as hydrogen 
bonding between neutral lipids, or indeed interactions between lipids across two apposed 
bilayers. 

5. Mixed solvents effects in ionic solutions 

So far, in most theoretical studies of interactions between charged macromolecular 
surfaces, the surrounding liquid solution was regarded as a homogeneous structureless 
dielectric medium within the so-called "primitive model". However, in recent 
experimental studies on osmotic pressure in solutions composed of two solvents (binary 
mixture), the osmotic pressure was found to be affected by the binary solvent 
composition [26]. It thus seems appropriate to generalize the PB approach of section 2 
by adding local solvent composition terms to the free energy. 

Our approach is to generalize the bulk free energy terms to include regular solution 
theory terms for the binary mixture, augmented by the non-electrostatic interactions 
between ions and the two solvents in order to account adequately for the preferential 
solvation effects [27] . In this generalized PB framework the mixture relative composition 
will create permeability inhomogeneity that will be incorporated into the electrostatic 
interactions. 

More specifically, our model consists of ions that are immersed in a binary solvent 
mixture confined between two planar charged interfaces. Note that the two surfaces are 
taken as homogeneous charge surfaces with negative surface charge a < 0. We do this 
to make contact with experiments on DNA that is also negatively charged [26] . Though 
the model is formulated on a mean-field level it upgrades the regular PB theory in two 
important aspects. First, the volume fractions of the two solvents, 0a and 0b = 1 — 0a, 
are allowed to vary spatially. Consequently, the dielectric permeability of the binary 
mixture is also a function of the spatial coordinates. In the following, we assume that 
the local dielectric response e{r) is a (linear) compositionally weighted average of the 
two permeabilities ^a and Eb- 




(27) 



§ A detailed discussion about the value of Xs is given in [21] 
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or, 

£(r) = £o-0(rK, (28) 

where we define cj) = (pB, £o = £a and = ea — £b- This linear interpolation assumption 
is not only commonly used but is also supported by experimental evidence [28]. Note 
that the incompressibility condition satisfies 0a + 0b = 1 , meaning that ionic volume 
fractions are neglected. 

The second important modification to the regular PB theory is the non-electrostatic 
short-range interactions between ions and solvents. Those are mostly pronounced at 
small distances and lead to a reduction in the osmotic pressure for macromolecular 
separations of the order 10-20 A. Furthermore, it leads to a depletion of one of the two 
solvents from the charged macromolecules (modeled here as planar interfaces), consistent 
with experimental results on the osmotic pressure of DNA solutions [26] . 

The model is based on the following decomposition of the free energy 

F = FpB + -Fmix + Fsol- 

The first term, Fp^ is the PB free energy of a 1:1 monovalent electrolyte as in with 
one important modification. Instead of a homogeneous permeability, e, representing a 
homogeneous solution, we will use a spatial-dependent dielectric function e{r) for the 
binary liquid mixture. The second term, -Fmix[0], accounts for the free energy of mixing 
given by regular solution theory: 



^ mix — ^3 I >-i ' 



'ln0+(l-0)ln(l-0) + x0(l-0^ 



The interaction parameter, x, is dimensionless (rescaled by k^T). As the system is 
in contact with a bulk reservoir, the relative composition has a chemical potential 
which is determined by the bulk composition 0^. For simplicity, we take the same 
molecular volume ~ for both A and B components. 

The third term, Fgoi, originates from the preferential non-electrostatic interaction of 
the ions with one of the two solvents. We assume that this preference can be described 
by a bilinear coupling between the ion densities, n±, and the relative solvent composition 
0. This is the lowest order term that accounts for these interactions. The preferential 
solvation energy, Fgoi, is then given by 

k^T / d^r (a+n+ + Q;_n_) (30) 



-Tsol 

where the dimensionless parameters a± describe the solvation preference of the ions, 
defined as the difference between the solute (free) energies dissolved in the A and B 
solvents. Finally, to all these bulk terms one must add a surface term, F,, [as in ([6])] 
describing the electrostatic interactions between charged solutes and confining charged 
interfaces. 

In thermodynamic equilibrium, the spatial profile of the various degrees of freedom 
characterizing the system is again obtained by deriving the appropriate Euler-Lagrange 
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Figure 5. Spatial profiles of (a) the solvent relative composition <f> and (b) the permeability 
e. The regular PB with homogeneous dielectric constant e — 77 (solid line) is compared 
with our modified PB for binary mixture with and without short-range interactions, 
X = (dashed line) and x ~ 1-5 (dotted line), respectively. Other parameters are: 
cr = -l/lOOA-2, nb = lO^'^M, EA = 80, Eb = 20 and </«b = 0.05. In all the cases, no 
preferential solvation is considered, a± — 0. 



(EL) equations via a variational principle. The EL equations are then reduced to four 
coupled differential equations for the four degrees of freedom, '?/'(r), n±(r) and 0(r). 
First we have the Poisson equation 

V ■ (—Vip) +e(n+-n) = (31a) 
V47r / 

then the Boltzmann distribution 

± + ln(n±a3) + a±(/) - /i± = (316) 
KbT 

and finally the EL equation for the density field 0(r) 

+ (a+n+ + a^n) - = (31c) 

At the charged interface, the electrostatic boundary condition stems from the variation 
of Fg with the difference that 6{(j)) has a surface value: Eg = Eq — e^(j)si so that the 
boundary condition becomes 

47re 

n ■ = a, (32) 

s Es 

Again the boundary condition states the electroneutrality condition of the system as 
can be shown by the integral form of Gauss' law. 

By solving the above set of equations, one can obtain the spatial profiles of the 
various degrees of freedom at thermodynamic equilibrium. For a general geometry. 
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Figure 6. The dependence of pressure on separation D for various ion-solvent interaction 
strengths a+. Other parameters are: a — — rtb = ea = 80, Eb — 4: 

and (pt = 0.09. 

these equations can be solved only numerically to obtain spatial profiles for and n±. 
The osmotic pressure can then be evaluated via an application of the first integral of 
(13 lap . (131 &|) . fl31cp .This pressure is a function of the inter-plate separation D and the 
experimentally controlled parameters a+, 0^ and nj,. 

For two identically charged planar surfaces, in the absence of preferential solvation, 
the numerical solutions of the EL equations show (see figure [5]) that the density profiles 
and, therefore, the osmotic pressure undergo only small modifications. By adding the 
preferential solvation term as quantified by «+, one observes a considerable correction 
to both the density profile, as well as the pressure (see figure |6]). Most notably is the 
reduction in the osmotic pressure at small separations (10-20 A) due to the coupling 
between ion density and solvent local composition. 

6. Polarizable ions in solution 

Ion-specific effects as manifested through the Hofmeister series have been recently 
associated with ionic polarizability, especially in the way they affect van der Wails 
electro dynamic interactions between ions and bounding interfaces [29]. However, this 
does not provide the full description, since the ionic polarizability also modifies the 
electrostatic interactions of ions with the surface charges. 

In what follows we generalize the PB theory in order to include also the contribution 
of the ionic polarizability to the overall electrostatic interactions. This inclusion leads 
to a new model that again supersedes the standard PB theory. Note, that in a more 
complete and consistent treatment, the polarizability should also be taken into account 
in the electrodynamic van der Wails interactions, in addition to what we propose here. 

We begin along the lines of section 5, while delimiting ourselves to the counterion 
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polarizability in the electrostatic part of the free energy. In order to discuss a manageable 
set of parameters, we discard other possible non-electrostatic terms in the total free 
energy and concentrate exclusively on the changes brought about in the ionic density 
profile and the interaction osmotic pressure. 

The free energy has the same form as FpB for counterion only case (|T2i) with one 
important difference that e = e{n) is now a function of the counterion density n 



F/A = / dz 



^p + emp + kBTn(\nn — 1) — fin 



(33) 



Sir 

The equilibrium equations read 

^ = ATien + — [e(nU'] = (34a) 
dtp dz 

6F ^ lde^^„ ^^^^ ^^^^^^ ^ ^ 
on Svr on 

The first equation is a generalization of the Poisson equation and the second is a 
generalization of the Boltzmann distribution. For the free energy (1331) we can use the 
general first integral of the system that gives the pressure in the form: 

P = -^(ein) + ^^n]tlj'^ + kBTn (35) 



^TT \ dn ^ 

The first term is an appropriately modified form of the Maxwell stress tensor while the 
second one is the standard van 't Hoff term. 

iFiom. the first integral we furthermore derive the following relation: 



Wn P) 
V + ^« 

where P = P/k^T is the rescaled pressure. Substituting this relation in (134 gp . we end 
up with a first-order ordinary differential equation for the ion density n: 

dn /27re2 n ^^^^ 



dz V ^bT df{n)/dn 



where 




/H=^HW , " a^w (38) 

is a function of the variable n only. Equation (1371) can be integrated explicitly either 
analytically or numerically, depending on the form of f{n). Note that in thermodynamic 
equilibrium, the total osmotic pressure P is a constant. 

The boundary condition for a constant surface charge is given by 

e{n,)i,'l = 47re|a| (39) 

where e{ns) and Ug are the surface values of the dielectric function and ion density, 

respectively. Using the pressure definition, we arrive at an algebraic equation for the 
surface ion density: 

-- 27re^(T^ / , , de(ns) \ 
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For a single surface the pressure vanishes P = (with analogy to two surfaces at infinite 
separation, D oo), and the basic equations simplify considerably. 

As a consistency check, we take the regular case, i.e. a single charged surface 
with a homogeneous dielectric constant. In this case the function f{n) takes the form 
f{n) = ^/eon, and (1371) reads: 
dn 



_ = -v/Svr/Bns (41) 

which gives the well known Gouy-Chapman result for the counter-ion profile close to a 
single charged plate (in absence of added salt) : 

^(^) = 1 \2 (42) 

2'kIb{z + Agc) 

where Aqc = l/(27r/B|cr|) is the Gouy-Chapman length, found by satisfying the boundary 
condition ( l39i) 

We assume that to lowest order the dielectric constant can be expanded as a function 
of the counterion concentration n as: 

e{n) = eo + (3n + 0{n^), (43) 

where Eq is the dielectric constant of the solvent and (3 = de / dn\o is a system parameter 
describing the molecular polarizability of the counterions in the dilute counterion limit. 
The derivative of /(n) to be used in ( 137|) reads: 

dl ^ 2PMn-P)+einMn)+Pn] ^^^^ 
{n - P)[e{n) + (3nY 

and (|371) can now be solved explicitly for n{z). 

In figure [7] we show the spatial profile of the counterion density away from single 
charged surface. No extra salt is added and the pressure is zero P = for a single plate 
as mentioned above. In figure [7^ we present the case where /5 < {i.e., the dielectric 
permeability in the regions of high ion density is smaller than in the pure solution). 
The ion density profile exhibits a somewhat "fiat" behavior, indicating some saturation 
in the vicinity of the surface, while at a distances of ~ 5 A and further away from the 
surface the density decays strongly to zero. In figure [7b the density profile is shown for 
the opposite case where /3 > 0. Namely, when the dielectric permeability in the regions 
of high ion density is larger than in the pure solution. The profile here does not show 
large deviations from the usual Gouy-Chapman behavior even close to the surface. 

Obviously, the effect of ionic polarizability strongly depends on the sign of the ionic 
polarizability, [3. It appears that the dependence of the dielectric constant on ionic 
density introduces effective interactions between the ions and the bounding surface. For 
/3 < 0, these additional interactions seem to be strongly repulsive and long ranged. 
They lead to a depletion of the ions in the vicinity of the surface. In the opposite case 
where for /3 > 0, the interactions are also repulsive but extremely weak and the electric 
double layer structure remains almost unperturbed by the ionic polarization. This sets 
a strong criterion for ion specificity because the ions can be differentiated according to 
the sign of their polarizability affecting their surface attraction. 
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Figure 7. Ion density and permeability profiles. The solid lines in both panels represent the 
case where the ions have no contribution to the permeability, [3 = 0, given by (j42|l . In (a) 
the dashed line shows the ion density where their averaged permeability is smaller than the 
solution dielectric constant, f3 = —30 M^^. The permeability of the second case is shown in 
the inset. In (b) the dashed line shows the ion density where their averaged permeability is 
smaller than the solution dielectric constant, P — 20 M^^. The permeability of the second 
case is shown in the inset. Other parameters in both (a) and (b) are: a = — 1/100A~^, 
and eo = 80. 



7. Conclusions 

We presented here several attempts to generalize the PB theory by including in the free 
energy additional terms to the standard electrostatic and ideal entropy of mixing. We 
also showed how these terms lead to modifications of the PB equation and its boundary 
conditions. More specifically, we have aimed to include additional interactions between 
dissolved ions, such as finite ion size and polarizability, as well as their solvation and 
interactions with bounding surfaces. All these endeavors are done within the mean-field 
approximation while neglecting charge density fluctuations and ion correlations. 

The approach presented here amends the free energy in specific ways giving rise to 
modified electric bilayer charge distribution and ensuing interactions between charged 
surfaces as mediated by ionic solutions. Some modification wrought by the non- 
electrostatic terms can be interpreted in retrospect as specific interactions of the ions 
with the bounding surfaces or between the ions themselves. 

To verify these models, results should be compared with appropriate experiments. 
In some cases, such comparisons are feasible and are qualitatively favorable [SHEZj. On 
the other hand, comparisons with extensive all-atom MC or MD simulations are not 
obvious since these simulations also contain multiple parameters that have no obvious 
analogue. It would, nevertheless, be valuable to link these approaches with appropriately 
designed experiments or even to check the predictions of these approaches in realistic 
systems. 
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